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Abstract 

Through Monte Carlo Simulation, the well-known majority-vote model has been 
studied with noise on directed random graphs. In order to characterize completely 
the observed order-disorder phase transition, the critical noise parameter qc, as well 
as the critical exponents ^jv and Xjv have been calculated as a function of 
the connectivity z of the random graph. 
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1 Introduction 



The majority- vote model (MVM) [T|2] . a nonequilibrium model defined by 
stochastic dynamics with local rules and with up-down symmetry, defined on 
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regular lattices shows a second-order phase transition with critical exponents 
f3,'y, u - which characterize the system in the vicinity of the phase transition - 
identical with those of equilibrium Ising model PPj . More general, it has 

been argued that the existence of up-down symmetry in two-state dynamic 
systems implies the same critical behavior (same universality class) of the 
equilibrium Ising model [71I31I8] . 

On the other hand, MVM on the complex networks exhibit different behav- 
ior, i.e., it belongs to different universality class Pll0flllll2j . Campos et al. 
investigated MVM on a small- world network [9], which was constructed using 
the square lattice (SL) by the rewiring procedure [13]. Campos et al. found 
that the critical exponents 7/z/ and l3/v are different from these of the Ising 
model [B] and depend on the rewiring probability. Pereira et al. [12] studied 
MVM on Erdos-Renyi's (ER) classical random graphs [131IT5] . and Lima et 
al. [To] also studied this model on random Voronoy-Delaunay [16] lattice with 
periodic boundary conditions. Very recently Lima [11] studied the MVM on 
directed Albert-Barabasi (AB) network [T7] and contrary to the Ising model 
[5j on these networks [TS], the order /disorder phase transition is observed in 
this system. The calculated l3/v and 7/1/ exponents are different from those for 
the Ising model [6] and depend on the value of connectivity z of AB network. 
The latter was observed also for undirected ER random graph [12] . 

In this paper we study the Majority-vote model with noise on directed Erdos- 
Renyi (ER) random graphs [T^fTS] . Through Monte Carlo (MC) simulations 
and standard finite-size scaling techniques we determine the critical exponents 
for several values of the connectivity z of the graph, in order to characterize 
the observed order-disorder phase transition. 



2 The Model 



The majority- vote model, on directed Erdos-Renyi (ER) random graphs, is de- 
fined [2|10|12|14II15II19] by a set of "voters" or spins variables = ±1, located 
on every site (or node) of a directed ER graph. The system dynamics is as 
follows: For each spin we determine the sign of the majority of its neighboring 
spins. With probability g, known as the noise parameter, the spin takes the 
opposite sign of the majority of its neighbors, otherwise it takes the same sign. 
It is important to emphasize that the noise parameter q plays the role of the 
temperature in equilibrium systems and measures the probability of aligning 
anti-parallel to the majority of neighbors. The probability Wi of a single spin 
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flip is given by 



1 - (1 - 2q)aiS[^ai+s 

■<5=1 



, with S{x) 



A, se X 7^ 
0, se X = 



where the summation is over all ki sites connected to spin o"j. It is worth to 
mention that this probability exhibits up-down symmetry, i.e., w{ai) = w{—ai) 
under the change of the spins in the neighborhood of (Tj. 

It is argued that on a square lattice there exists two phases for sufficiently 
small q, with a formation of an island of up spins on a sea of down spins. The 
size of this island follows a birth-and-death process in which the death rate is 
larger than the birth rate, thus preventing the growth of the island and keeping 
the down spin phase stable. By symmetry, there must exist another phase 
with spins up. However, if the up-down symmetry of the spin flip probability 
is broken, no phase transition is expected, as in the Ising model [1]. In two 
dimensions this model has a ferromagnetic stationary phase for < g < gc 
undergoing a second-order phase transition to a paramagnetic phase at qc 
{qc = 0.075 for a square lattice). The static critical behavior is Onsager-like 
[2][3ll8] . and according to the argument of Grinstein et al. its dynamic 
critical behavior is the same as model Ising. 

Indeed, the majority-vote is a particular case of a general class of polling 
models [201I2T] composed by interacting two-state opinion agents. The main 
result of the most opinion formation models is that the dynamical rule leads to 
an opinion polarization of the whole population along one of the two competing 
opinions [20ll21f22j . Recently, the contrarian concept - agents which have the 
opinion opposite to that of the majority - was introduced to take into account 
some peculiar behavior of the agents [23] • With a small fraction of contrarians, 
the system no longer leads to total polarization: the minority opinion does 
not disappear from the population. Above a critical fraction of contrarians, 
polarization does not occur, and the two possible opinions are shared each by 
half the agents [23ll2i] . 

Depending on the nature of the interactions, networks can be directed or 
undirected [25]. In directed networks, the interaction between any two nodes 
has a well-defined direction, which represents, for example, the direction of 
material fiow from a substrate to a product in a metabolic reaction, or the 
direction of information fiow from a transcription factor to the gene that it 
regulates. In undirected networks, the links do not have an assigned direction. 
For example, in protein interaction networks a link represents a mutual binding 
relationship: if protein A binds to protein B, then protein B also binds to 
protein A. In other words, in undirected networks, if a node A is linked to B, 
then B must also be linked to A, while in directed ones, if a node A is linked 
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to B, B may not be linked to A, but to some node else instead. In the present 
work, we only consider directed graphs. 



Starting with disconnected nodes, the Erdos-Renyi network are generated 
by connecting couples of randomly selected nodes with a probability < p < 1 
for every pair, prohibiting multiple connections (i.e., couples of nodes con- 
nected by more than one link), until the number of edges equals z (connec- 
tivity or number of links of a node). In the limit — > oo, it is found that the 
tail (high z region) of the degree distribution P{z) decreases exponentially, 
which indicates that nodes that significantly deviate from the average are ex- 
tremely rare. The clustering coefficient - which is the probability that two 
neighbors of the same node are also mutual neighbors [TH|E5] - is independent 
of a node's degree, so C{z) appears as a horizontal line if plotted as a function 
of z. The mean path length is proportional to the logarithm of the network 
size, / ~ logN [15], which indicates that it is characterized by the small- world 
property. 



To study the critical behavior of the model we consider the magnetization M , 
the susceptibility x and the Binder's fourth-order cumulant, which are defined 
by 



N 



UN[q) = 1 



3((m2),), 



(2) 
(3) 
(4) 



where N is the number of vertices of the random graph with fixed z, (. . .)t 
denotes time averages taken in the stationary regime, and (. . .)s stands for 
sample averages (20 samples). 

These quantities are functions of the noise parameter q and, in the critical 
region, obey the following finite-size scaling relations 



M = iV-^/^/^(x)[l + ...], 



(5) 
(6) 



dU 
dq 



iVV-/^(x)[l + 



(7) 
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the brackets [1 + ...] indicate corrections-to-scaling terms, i/, (3 and 7 are the 
usual critical exponents, fi{x) are the finite size scaling functions with 



X 



{q - q,)N^/^ (8) 



being the scaling variable. Therefore, from the size dependence of M^v and x 
we obtained the exponents (3/v and 7/z/, respectively. The maximum value of 
susceptibility x ^Iso scales as N^^'^. Moreover, the value of q for which x has 
a maximum, gj™="' = qc{N), is expected to scale with the system size as 

gc(iV) = gc + &iV-v-, (9) 



were the constant b is close to unity. Therefore, the equations [7] and [9] are used 
to determine the exponent l/iy, as well as to check the values of qc obtained 
from the analysis of the Binder's cumulant (Eq. H]). Finally, we have also 
examined if the calculated exponents satisfy the hyperscaling hypothesis 

2/3/z/ + 7/z/ = Deff (10) 



in order to get the effective dimensionality, -Dgg, for various values of connec- 
tivity z. 

We have performed Monte Carlo simulation on directed ER random graphs 
with various values of connectivity z. For a given z, we used systems of size 
N = 250, 500, 1000, 2000, 4000, 8000 and 16000. We waited 10000 Monte Carlo 
steps (MCS) to make the system reach the steady state, and the time averages 
were estimated from the next 10000 MCS. In our simulations, one MCS is 
accomplished after all the N spins are updated. For all sets of parameters, we 
have generated 20 distinct networks, and have simulated 20 independent runs 
for each distinct network. 



3 Results and Discussion 



In Fig. [T] we show the dependence of the magnetization and the suscep- 
tibility X oil the noise parameter, obtained from simulations on directed ER 
random graphs with N = 16000 sites and several values of connectivity z. 
In Fig. 1(a) each curve for Mjy, for a given value of and z, suggests that 
there is a phase transition from an ordered state to a disordered state. The 
phase transition occurs at a value of the critical noise parameter q^ which is 
an increasing function the connectivity z of the directed ER random graphs. 
In Fig. 1(b) we show the corresponding behavior of the susceptibility x- The 
value of q where x has a maximum is here identified as qdN) (Eq. [2]). 
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Fig. 1. (a) Magnetization and (b) susceptibility as a function of the noise parameter 
q, for = 16000 sites. From left to right we have z = 4, 6, 8, 10, 20, 50, and 100. 




Fig. 2. Binder's fourth-order cumulant as a function of q for several values of network 
sizes N {N = 250, 500, 1000, 2000, 4000, 8000 and 16000), and for two different 
values of the connectivity z: z = A (on the left-hand side) and z = 50 (on the 
right-hand side). 

To determine estimates for the critical point gc, we calculate the Binder fourth- 
order magnetization cumulant U at different values of the noise q and several 
network sizes A^. Finite size scaling predicts that for sufficiently large systems, 
these curves should have a unique intersection point U* [27]. The value of q 
where this crossing occurs is the value of the critical noise qc which is not 
biased by any assumptions about critical exponents, since by construction, 
the Binder cumulant presents zero anomalous dimension, therefore it respects 
the correct critical behavior of the system near PTPS] . 

In Fig. [2] we plot the Binder's fourth-order cumulant for different values of 
and two different values of z {z = 4 and z = 50). The critical noise parameter 
qc, for a given value of z, is estimated as the point where the curves for 
different system sizes N intercept each other. As we can noticed, there exists 
some dependence of the critical noise qc with the connectivity z, i.e., when the 
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(a) 



(b) 



Fig. 3. (a) The intersection point U* versus the connectivity z; (b) Phase diagram: 
the critical values of the noise parameter qc as a function of the connectivity z. 

connectivity z increases, the critical noise also increases: for 2; = 4, gc ~ 0.18 
(left-hand side of accept Fig. [2]), and for z = 50, qc ~ 0.412 (right-hand side 
of Fig. El). This same procedure is also made for several values of connectivity 
z, in this way the Binder's cumulant crossing point U* as a function of the 



connectivity z could be obtained (Fig. 3(a)), as well as the phase diagram is 
built, i.e., the dependence of the critical noise qc with the node connectivity 
z, what is shown in Figure 3(b), Based on the fact that the crossing point 



of the cumulant (for different system sizes) gives the transition point and the 
value of the cumulant U* at the transition point indicates the universality class 
of the transition [27]. So, from Fig. 3(a), which shows that the intersection 
point U* of the cumulant decreases when the connectivity z increases, it can 
be argued the present model does not fall into the same universality class of 
the equilibrium Ising model, since U* = 2/3 well inside the ordered phase 
and U —>■ well within the disordered phase. The Fig 3(b) indicates the 
increase of the critical noise qc with the increase of the node degree z, and 
moreover, for higher values of z, the critical noise qc approaches to that one 
obtained {qc ~ 0.5) from mean-field theory when z — > cxo [12]. In addition, 
according to the mean- field theory [T2|, in the limit the connectivity z goes to 
infinity, the magnetization approaches also to the value 0.5, what one could 
infer from Fig [H It is important to emphasize that these results obtained 
from our Monte Carlo simulations are in a good agreement with a previous 
work of the majority- vote model with noise on random graphs [12], in which 
simulations and mean-field analysis were performed to characterize the system 
transition. 



Figure 4(a) and Figure 4(b) show, at g = q^ the dependence of the mag- 



netization and the susceptibility with the system size A^, respectively, when 
different values of the connectivity z are considered. As we can see, the ob- 
tained straight lines, whose slopes correspond to the exponents ratio P/u (Fig. 



4(a)) and 7/1/ (Fig. 4(b)), confirm the scaling equations, given by Eqs. Oand 
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In N 



(a) 



In N 



(b) 



Fig. 4. (a) In M{qc) versus In N and (b) In x{Qc) versus In A^. In both figures, from 
top to bottom, z = 4, 6, 8, 10, 20, 50 and 100. 

[6l Moreover, while there is a shght tendency for the exponent P/u to increase 
with z, the opposite occurs to the exponent 7/z/ that decreases with z. 



2 




10 



In N 



Fig. 5. Log- log plot of the susceptibility at its maximum Xn 
(circles) versus N, for connectivity 2: = 8. 



(squares) and q = Qc 



In Fig. 0, we present the dependence, at g = qc, of the susceptibility x of 
its maximum value Xmax with the system size N. The exponent ratio 7/1/ is 
obtained from the slopes of the straight lines. For almost all the values of z, 
the exponents 7/1/ of the two estimates agree within error bars (see Tabled]). 



As already observed in Fig. 4(b) , an increased z means a slight tendency to 



decrease the exponent ratio 7/z/. 



z 










Xjv 




4 


0.175(4) 


0.230(5) 


0.530(6) 


0.516(2) 


0.545(26) 


0.990(7) 


6 


0.238(3) 


0.243(4) 


0.509(4) 


0.511(2) 


0.488(16) 


0.995(5) 




0.z74(3j 


c\ 000/ a\ 

0.238(4) 


0.512(4) 


0.504(2) 


0.548(14) 


0.988(5) 


10 


0.299(2) 


0.268(4) 


0.473(5) 


0.495(1) 


0.487(10) 


1.009(6) 


20 


0.359(2) 


0.280(4) 


0.451(4) 


0.485(2) 


0.510(10) 


1.011(5) 


50 


0.412(2) 


0.282(3) 


0.441(2) 


0.466(5) 


0.484(11) 


1.005(3) 


100 


0.438(2) 


0.286(2) 


0.428(3) 


0.440(8) 


0.520(19) 


1.000(3) 



Table 1 

The critical noise (/c? cind. the critical exponents, for directed ER. random graphs with 
connectivity z. Error bars are statistical only. 

The Table [H summarizes the values of Qc, the three critical exponents jS/u, 
7/z/, l/u and the effective dimensionality of the system. For all values of z, 
= 1, which has been obtained from the hyperscaling hypothesis (Eq.lO), 
therefore when /3/z/ increases, 7/1/ decreases at qc, thus providing Defr = 1 
(along with errors). 

It is worth to mention that while our results for the majority-vote model on 
a directed ER random graph are in good agreement with the model on other 
complex networks, they are different from those found for the majority-vote 
model and the traditional Ising model on a square lattice Moreover, the 
critical exponents and the hyperscaling relation obtained from our simula- 
tions corroborate a proposition that the majority-vote models defined on a 
regular lattice, on small- world networks [9], on Voronoy-Delaunay lattice [TO] , 
on Barabasi- Albert networks [TT][29] . and on undirected ETdos-KenyVs random 
graphs [12\, belong to different universality classes. 



4 Conclusion 

Using Monte Carlo simulations, we found a second-order phase transition of 
the majority-vote model with noise on directed ER random graphs. The ob- 
served phase transition, which occurs with connectivity z > 1, was completely 
characterized through the phase diagram and the critical exponents. Although 
the obtained exponents are different from those for the same model on other 
topologies [2|9|10fllll29f30j . they suggest that the majority- vote models de- 
fined on a regular lattice, on small-world networks, and on Erdos-Renyi's 
random graphs, belong to different universality classes. In fact, it was found 
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that the critical exponents depend on the long-range interactions (short cuts), 
i.e., they depend on the node degree 2;, in such way that results in a larger 
robustness of the system against the noise: the larger the connectivity 2, the 
greater the critical noise Qc- The found values U* of the Binder's fourth-order 
cumulant at the transition point for different values of the connectivity z 
were different from those obtained for the isotropic majority-vote model on 
a regular square lattice and, moreover, for the square Ising model, what cor- 
roborates the conjecture that the majority-vote model on complex network 
does not belong to the same universality class of the equilibrium Ising model. 
Finally, we should remark that the importance of our results is not only in 
presenting a set of numbers, but it is in characterizing completely the found 
phase transition, i.e., to determine the whole set of the critical exponents of 
the order- disorder phase transition appearing in the system. Moreover, the full 
set of critical exponents for several values of the connectivity z, which is really 
different from the exact values of the Ising model and majority-vote model on 
a regular lattice, shows clearly the influence of the small-word effects in the 
majority-vote model, as well as it points out the dependence of the critical 
exponents with the connectivity z. Therefore, the well-known effects of the 
shortest path length in systems with small-world topology on the information 
propagation are the responsible for making more effective the influence of the 
neighborhood, which now is not constrained by the geographical distance be- 
tween the individuals, as observed in regular lattices. Of course, we are aware 
that random graphs are seldom very good models for the type of networks one 
finds in nature [T7I3T]- While they have a low shortest path length, they do 
not have another important property of most observed networks: clustering. 
This is obviously not the case in our studied topology, where all vertex pairs 
are independently connected with the same probability. 
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